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We introduce a definition of a "localization width" whose logarithm is given by the entropy of the 
distribution of particle component amplitudes in the Lyapunov vector. Different types of localization 
widths are observed, for example, a minimum localization width where the components of only two 
particles are dominant. We can distinguish a derealization associated with a random distribution of 
particle contributions, a derealization associated with a uniform distribution and a derealization 
associated with a wave-like structure in the Lyapunov vector. Using the localization width we show 
that in quasi-one-dimensional systems of many hard disks there are two kinds of dependence of the 
localization width on the Lyapunov exponent index for the larger exponents: one is exponential, and 
the other is linear. Differences, due to these kinds of localizations also appear in the shapes of the 
localized peaks of the Lyapunov vectors, the Lyapunov spectra and the angle between the spatial 
and momentum parts of the Lyapunov vectors. We show that the Krylov relation for the largest 
Lyapunov exponent A ~ — plnp as a function of the density p is satisfied (apart from a factor) in 
the same density region as the linear dependence of the localization widths is observed. It is also 
shown that there are asymmetries in the spatial and momentum parts of the Lyapunov vectors, as 
well as in their x and i/-components. 

Pacs numbers: 05.45. Jn, 05.45.Pq, 02.70.Ns, 05.20.Jj 
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I. INTRODUCTION 

The dynamical instability of a system is essentially 
connected to the unpredictability of that system. It is 
described by the time-evolution of the difference of two 
phase space vectors describing nearby trajectories, the so 
called the Lyapunov vector. If the amplitude of the Lya- 
punov vector increases (decreases) rapidly in time, then 
the dynamics is unstable (stable) in the direction of the 
Lyapunov vector. An unstable orbit implies that a part 
of the dynamics is unpredictable and a statistical treat- 
ment may be required. The Lyapunov exponent, defined 
as the rate of exponential divergence or contraction of 
the amplitude of the Lyapunov vector with time, is the 
most frequently used indicator of the dynamical insta- 
bility, and a system with a nonzero positive Lyapunov 
exponent is called chaotic. The Lyapunov exponent is 
connected to the transport coefficients, such as the con- 
ductance and viscosity [lj . In general, for each Lyapunov 
vector there is an individual Lyapunov exponent that is 
different in magnitude, and this leads to the concept of a 
sorted set of Lyapunov exponents, the so called Lyapunov 
spectrum. 

Although the significance of the amplitudes of the Lya- 
punov vectors has been emphasized in the discussion of 
the dynamical instability, the Lyapunov vector itself in- 
cluding the information about its anqle can play an im- 
portant role in chaotic dynamics [J, 0. For example, 
the Lyapunov vector was used to characterize a high- 
dimensional chaotic attractor and the clustered mo- 
tion in symplectic coupled map systems Q. One may 
also mention that the wavelikc structure of the Lyapunov 
vectors associated with the stepwise structure of the Lya- 
punov spectra in the small (in absolute value) Lyapunov 



exponent region has been reported in many-hard-disk 
systems 0,00 E3 

The localization of the Lyapunov vector for high- 
dimensional chaos, which we call the "Lyapunov localiza- 
tion" for convenience in this paper, is one of the chaotic 
behaviors which involves information on all the individ- 
ual components of the Lyapunov vectors. The Lyapunov 
localization appears as the spatial localization of the 
largest components of the Lyapunov vector and is gener- 
ally associated with the large Lyapunov exponents. Such 
a phenomenon is reported in the Kuramoto-Sivashinsky 
model , in coupled map lattice systems 0, , in a 
random matrix model |12| . in a high-dimensional sym- 
plectic map system |13| and in many-hard-disk systems 
|l4l Il5j . However it was not so clear whether the Lya- 
punov localization has its origin in a randomness pro- 
duced by the chaotic dynamics 0] (like Anderson local- 
ization [13) or it comes simply from a short range prop- 
erty of particle interactions. Moreover physical meaning 
and significances of this phenomenon has not been well 
understood. One of the few suggestions about the im- 
portance of the Lyapunov localization is that it may be 
related to the existence of the thermodynamic limit for 
the largest Lyapunov exponent ^4j. If a spatial local- 
ized region of the Lyapunov vector corresponding to the 
largest Lyapunov exponent is independent of the particle 
number N in the thermodynamic limit, then the largest 
Lyapunov exponent can be ./V-independent. This implies 
that the Kolmogorov-Sinai entropy, which is equal to the 
sum of all the positive Lyapunov exponents in closed sys- 
tems, is an extensive quantity like the thermodynamic 
entropy. Besides, this gives some supporting evidence 
to the existence of the thermodynamic limit of the Lya- 
punov spectrum, which has been a subject of study in 
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many-particle chaotic systems 0, 0, 0, Hjj . 

The principal aim of this paper is to use the Lyapunov 
localization as an indicator of the chaotic properties of 
many-particle systems. For quantitative considerations 
of the Lyapunov localization we use an entropy-like quan- 
tity (or information dimension) of a distribution function 
of the amplitudes of the Lyapunov vector components 
for each particle [Til fl3| . We introduce the " localization 
width" as the quantity whose logarithm is given by this 
entropy- like quantity. The value of the localization width 
is in the range [1, A], and indicates the number of par- 
ticles contributing to the localized part of the Lyapunov 
vector components. Using the localization width we can 
also distinguish different types of delocalized behaviors 
of the Lyapunov vectors such as a derealization associ- 
ated with a random distribution of particle component 
amplitudes of the Lyapunov vector, a derealization as- 
sociated with a uniform distribution and a derealization 
associated with a wave-like structure in the Lyapunov 
vector. As a concrete system to consider for the study 
of Lyapunov localization, we use a quasi-one-dimensional 
system consisting of many hard-disks, in which the sys- 
tem shape is so narrow as to exclude exchange of parti- 
cle positions. In this system the minimum value of the 
localization width is given by 2, because particle interac- 
tions are given by collisions between two particles. In the 
quasi-one-dimensional systems each particle can interact 
only with its nearest neighbor particles, so the numer- 
ical calculation is less time-consuming than for a fully 
two-dimensional hard-disk system in which each parti- 
cle can collide with any other particle. The quasi-one- 
dimensional system also has an advantage in that the 
roles of the coordinate directions are strongly separated. 
In our system the narrow direction (the y-direction) of 
the rectangle is very much different from the longer or- 
thogonal direction (the x-direction) . Another advantage 
of the use of the quasi-one-dimensional system is that 
there is a wider region of stepwise structure in the Lya- 
punov spectrum compared to a square two-dimensional 
system |2l|. This is a noticeable point because in this 
paper we show that the localization width is an indi- 
cator not only of the Lyapunov localization but also of 
the stepwise structure of the Lyapunov spectrum. Cal- 
culating the localization width in quasi-one-dimensional 
systems we show that there are two kinds of dependence 
of the Lyapunov index on the localization widths: an 
exponential dependence [£ T>] and a linear dependence 
[CD]. (Here we define the Lyapunov spectrum as the set 
{AW, \( 2 \ • • • , A( 4JV )} of the Lyapunov exponents satis- 
fying the condition A^ > > ■ ■ ■ > X i4N) and in- 
troduce the the Lyapunov index as the ordering num- 
ber of the exponents in the Lyapunov spectrum.) The 
exponential dependence [£D] of the localization width 
appears at any particle density, and shows a tail in the 
spatial localized shape of the Lyapunov vector. On the 
other hand the linear dependence [CD] of the localization 
widths as a function of Lyapunov index appears only in 
cases of low particle density, and is characterized by a 



sharp rectangular localized shape of the Lyapunov vec- 
tor in space. Next, we consider the effects of the these 
two kinds of Lyapunov localizations on other quantities, 
such as the shape of the Lyapunov spectra, the angle be- 
tween the spatial part and the momentum part of the 
Lyapunov vector, and the amplitudes of the spatial part 
of the Lyapunov vector, etc. In particular, it is shown 
that apart from a prefactor the Krylov relation for the 
largest Lyapunov exponent A ~ —plrtp as a function of 
the particle density p is satisfied in the same density re- 
gion in which the linear dependence [CD] of the local- 
ization widths appears. We also demonstrate a relation 
between localized regions of the Lyapunov vectors and 
the positions of colliding particles, suggesting that the 
Lyapunov localization comes from the short range inter- 
action of the particles. We compare some of our results 
in quasi-one-dimensional systems with the square system 
with the same area, and show that the two kinds of lo- 
calizations also appear in the square case as well. 

As the second aim of this paper we investigate how the 
particle density and system shape affect the Lyapunov 
vector components. We show that the spatial part and 
the momentum part of the Lyapunov vectors are in al- 
most the same direction at high density, whereas they 
are rather close to orthogonal in low density cases, espe- 
cially in the region of the exponential dependence [£D] 
of the localization widths as a function of Lyapunov in- 
dex. The amplitudes of the spatial part of the Lyapunov 
vectors are larger (smaller) than the corresponding mo- 
mentum part in low density cases (high density cases). 
We also demonstrate that gaps appear in the amplitudes 
of the ^-components and the y-components of the spatial 
part and the momentum part of the Lyapunov vectors, 
because of the difference of the roles of the directions in 
the quasi-one-dimensional system. In these amplitudes of 
the Lyapunov vector components we can also see effects 
of the stepwise structure of the Lyapunov spectra. 

The outline of this paper is as follows. In Sect, [n] 
we compare some quantities that characterize the Lya- 
punov localization, and discuss the relative merits of the 
localization width comparing it with the other quanti- 
ties. The relation between the localized region of Lya- 
punov vectors and the position of colliding particles is 
demonstrated. In Sect. IIIII we show that there are two 
kinds of Lyapunov localizations. These two Lyapunov 
localizations are distinguished by their Lyapunov index 
dependences and the shapes of the localized peaks of the 
Lyapunov vectors. In Sect. II VI we investigate the effects 
of the two kinds of Lyapunov localizations on the shape of 
the Lyapunov spectra, and angles and amplitudes of the 
Lyapunov vectors components, etc. The Lyapunov in- 
dex dependence of the x and y-components of the spatial 
part and the momentum part of the Lyapunov vectors 
are shown. In Sect. [V]we investigate the density depen- 
dences of quantities associated with the largest Lyapunov 
exponent such as the largest Lyapunov exponent itself, 
the angle of the Lyapunov vectors and the localization 
width, etc., and specify the density region in which the 
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linear dependence [CD] of the localization width on Lya- 
punov index appears. The region of Lyapunov indices, in 
which the linear dependence [CD] of the Lyapunov widths 
appears, is connected to the region in which the Krylov 
relation for the density dependence of the largest Lya- 
punov exponent is satisfied. Sect. IVII is our conclusions 
and further remarks . In Appendix^Jwe give a derivation 
of an inequality for the localization width. In Appendix 
[H]we compare some results for the quasi-one-dimensional 
systems in the text with those of corresponding square 
systems. 



II. LOCALIZATION WIDTH OF LYAPUNOV 
VECTORS 

We consider the quasi-one-dimensional system which 
consists of N hard-disks with radius R. The shape of 
the system is a two-dimensional rectangle satisfying the 
condition 



RNV3 < L x and 2R < L v < 4R, 



(1) 



with width L x and height L y . Condition JQ| means that 
the system is so narrow that the positions of particles 
cannot be exchanged. The systems satisfying this condi- 
tion are referred to as quasi-one-dimensional systems in 
this paper. The schematic illustration of such a system 
is given in Fig. ^ The quasi-one-dimensional system 
was used to discuss the stepwise structure of the Lya- 
punov spectra and the associated wave-like structure of 
the Lyapunov vectors in Ref. [2l| . 



FIG. 1: Quasi-one-dimensional system. The system shape is 
so narrow and the exchange of particle positions are excluded. 



In this paper we consider the case where the parameter 
values are given by; the radius of particles is R — 1, the 
mass of particles is M = 1 , and the total energy of system 
is E = N. The system size is given by L y = 2R(l + 10~ 6 ) 
and L x = NL y (\ + d) satisfying the condition with 
a constant d, (except in Appendix [5] where we consider 
square systems.) Roughly speaking the factor 1 + d is 
the averaged ratio of the distances that each particle can 
move in the x-direction and y-direction. The constant d 
is connected to the density 



many-hard-disk systems the free flight part of the dy- 
namics is integrable, so the actual numerical calcula- 
tions of phase space dynamics and tangent space dynam- 
ics are simply described by multiplications of the time- 
evolutional matrices corresponding to free flight and col- 
lision. Second, in the quasi-one-dimensional system de- 
fined by Eq. Q each particle can only collide with its two 
nearest-neighbor particles, so we do not need to search 
every particle pair to find its colliding particles in numer- 
ical calculations. These points lead to a faster numer- 
ical calculation of the dynamical properties (especially 
Lyapunov spectra and Lyapunov vectors) than for other 
many-particle systems, such as fully two-dimensional 
particle systems (in which any pair of particles can col- 
lide), or with particles with soft core interactions. As 
another advantage of the quasi-one-dimensional system 
is that the roles of the x and y-directions are separated, 
so we can investigate how such a separated role for the di- 
rections effects for example, the Lyapunov exponents and 
Lyapunov vectors. In the quasi-one-dimensional system 
the information of the particle positions can be roughly 
replaced by the particle indices, so we can discuss the dy- 
namics of quantities by using three-dimensional graphs as 
functions of collision number and particle index, which is 
much simpler than the fully two-dimensional system that 
requires four-dimensional graphs as functions of the time, 
the x and y-components of particle positions. Finally in 
quasi-one-dimensional many-hard-disks, we get a wider 
stepwise region of the Lyapunov spectra rather than in 
the fully two-dimensional square systems. This was an 
essential point in the study of the stepwise structure of 
Lyapunov spectra and the associated wave-like structure 
of the Lyapunov vectors in systems with small numbers 
of particle in Ref. |2lj . 

We introduce the Lyapunov vector as ST^ n '(t) = 



(ST^i^JT^it), ■■■ , 5T { ">(t)) corresponding to the n- 

th Lyapunov exponents at time t. Here, STj (t) 
is the Lyapunov vector component corresponding to the 
j-th particle and the Lyapunov exponent A*™) at time t. 
We define the normalized amplitude 7^ (t) of the Lya- 
punov vector components ST^ (t) by 



(«)/ 



so that the conditions 



Eti *n nj (*) 



(«)/ 



(3) 



NttR 2 

L X Ly 



ttR 



(l + d)Ll 



(2) 



as d = irR 2 /(pLl) - 1 H. 

The quasi-one-dimensional many-hard-disk system has 
many advantages for numerical investigations of the dy- 
namical properties of many-particle systems. First, in 



N 

0< 7l (n) (t)<l 



(4) 



(5) 



are automatically satisfied. Fig. EJa) and (b) are two 
typical behaviors of the quantity 7^(i) corresponding 
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to largest Lyapunov exponent as functions of the colli- 
sion number nt and the particle index j in the case of (a) 
d = 1.5 (or density p = 0.314 • • • ) and (b) d = 10 5 (or 
density p — 0.00000785 ■ • •). Here we take the particle 
index j so that the ^-component q X j (t) of the j-th parti- 
cle satisfies the condition q x \{t) < q X 2 (*)<■•• < q X N(t), 
and the data is taken every two collisions. (It should be 
noted that the collision number nt is related to the time 
t approximately by multiplying by the mean free time, 
and in the quasi-one-dimensional system the particle in- 
dex j has a similar meaning to the x-component of the 
particle position.) These graphs show clearly that a non- 
zero value of the quantity (t) is localized in a small 
spatial region involving a few particles. Localized peak 
positions stay in almost the same position over several 
tens of collisions, then seems to hop to another position. 
Detail of the localized behaviors of the quantity 7 j (t), 
especially the differences between Figs. Ufa) and (b), will 
be discussed in the following sections and are the main 
purpose of this paper. 



be replaced by the limit of the amplitude of the spatial 
coordinate.) This quantity characterizes the Lya- 
punov localization as its spatial exponential decay rate. 
However this quantity may not be suitable to charac- 
terize the localization like that in Fig. Hfb) which does 
not have an exponential decay. Besides, this quantity 
requires the thermodynamic limit N — > oo. The numeri- 
cal calculation of the Lyapunov spectrum and Lyapunov 
vectors for many-particle systems is very time-consuming 
and has so far only been reported for systems of about 
1000 particles ^5|, so it may be rather difficult to esti- 
mate the quantity defined by Eq. JfjJ) numerically. 

Another method to characterize the Lyapunov local- 
ization was proposed in Ref. 0]. I* 1 this method we 
first introduce a parameter threshold 8 <E (0, 1) and de- 
fine the quantity C@ as 



Nx 



Cq = min ■ 



i=i 



~(n) 

r 3 



(*) 



(7) 



FIG. 2: Localized behaviors of the amplitudes 7!! of the 
normalized Lyapunov vector components of the j-th parti- 
cle corresponding to the largest Lyapunov exponent A' 1 ^ as 
functions of the collision number nt and the particle index j 
in a quasi-one-dimensional system of TV = 50. (a) The high 
density case of d = 1.5. (b) The low density case of d — 10°. 



To obtain our numerical results we use the algorithm 
developed by Benettin et al. |2^| and Shimada et al. [24[ 
and so on. This algorithm is characterized by intermit- 
tent rescaling and renormalization of Lyapunov vectors, 
which can be done after each particle collision in a many- 
hard-disk system. Other articles such as Refs. |25ll26. 27] 
should be referred to for more details of this algorithm 
and the Lyapunov vector dynamics of many-hard-disk 
systems. 

Now we discuss methods to characterize the strength 
of the localization of the Lyapunov vectors like those in 
Fig. |2ta) and (b) in a quantitative sense. By analogy to 
the localization length used in Anderson localization 
it may be suggested that the strength of the Lyapunov 
localization can be characterized by a localization length 
fi(") defined by 



with an integer Nx and the sorted set {7 1 ™' ) (i), 72™ "*(£), 
■■■,lP(t)} = {7{ n) (t),7^ n) (t), ■••,7&°(i)} satisfying 
the condition j[ n) (t) > j { 2 n \t) > ••• > 7 ^ } (t). This 

(n) 

quantity Cq is a measure of the number of particle com- 
ponents actively contributing to a localized part of the 
Lyapunov vector for the n-th Lyapunov exponent X^ n '. 

(n) 

This quantity C e does not require the thermodynamic 
limit, and can be suitable for both types of the Lyapunov 
localizations shown in Fig. [21 However this quantity is 
a function of an artificial threshold O, which cannot be 
determined by the dynamics itself. 

In this paper we discuss the Lyapunov localization by 
using an entropy-like quantity for the amplitude distribu- 
tion 7j"'(i) of the normalized Lyapunov vector elements, 
namely the entropy-like quantity defined by 



S (n) = 



N 

■E 



(n) /,s q (n) 

7} TO In 7} 



noting that the quantity 7 j n (t) satisfies the conditions 
(@J and ||SJ| so it can be regarded as a kind of proba- 
bility distribution function. Using this quantity we 
introduce the quantity as 



(*) 



(8) 



= lim lim 



ln 7 j n) (t) 



(6) 



yy(«) 



exp 



{S (Tl) }, 



(9) 



as the number of particles N goes to infinity, namely in 
the thermodynamic limit. Here we use the bracket (X) 
to signify the long time average of the time-dependent 
quantity X. (Note that in the definition © of the local- 
ization length f2( n ) we used the fact that the system is 
quasi-one-dimensional. In more general cases like a fully 
two-dimensional system the limit j — *■ 00 in Eq. JBJ) must 



which we call the n-th " localization width" corresponding 
to the n-th Lyapunov exponent A*-™) . This localization 
width has already been used to discuss the Lyapunov 
localization in a coupled map lattice model pd| and a 
high-dimensional symplectic map system |l3|. 

To understand the physical meaning of the localization 
width y\)( n ) defined by Eq. © it is useful to discuss some 
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of its properties. The first property of the localization 
width is the inequality 



1 < W {n) < N. 



(10) 



The derivation of the inequality i|l(J|) is given in Appendix 
lAl The second property of the localization width W( n ) is 
that the equality W^") = 1 is satisfied only when one of 
the quantities 7^ (t), j = 1, 2, ■ • • N is equal to 1, namely 
in the most localized case, and the equality = N 

is satisfied only when all of these quantities are equal to 
each other, namely j[ n \t) — 73 (t) = ■ ■ ■ = 7^ (t) = 
1/N, namely in the most delocalized case (See Appendix 
about this point.). These properties suggest that the 
localization width VVW also indicates the number of par- 
ticles contributing to the localized part of Lyapunov vec- 
tor, which is analogous to the quantity NC^ derived 
from Eq. J7J. The third property of the localization 
width comes from the symplectic structure of the Hamil- 
tonian dynamics. To discuss this property we note the 
conjugate relation [IJ 



ST 



(4JV-n+l) _ 



Sq 
Sp 



3 

(4JV-?i+l) 



X (") 

Sp) 

X (™) 



(ii) 



for the Lyapunov vector component ST^ 4N n+1 ^ of the j- 
th particle corresponding to the (4A^— n+l)-th Lyapunov 
exponent 



^(AN-n+1) _ _^(n) 



(12) 



n = 1, 2, • • • , 2N, where Sq^ and Sp^ are the spatial 
part and the momentum part of the Lyapunov vector 
ST^ of the j-th particle, respectively, and £M is a con- 
stant depending on the Lyapunov index only. (As a re- 
mark, if we use the Benettin algorithm to calculate the 
Lyapunov spectra, which is characterized by the intermit- 
tent normalization (as well as the rescaling) of Lyapunov 
vectors, the factor ^ n > in Eq. I|ll|) is given by +1 or —1 
depending on the Lyapunov index n.) From Eqs. © and 
(fTT|> we derive 
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(4JV— n+l) 



(*) 



(13) 



at any time t, where j is the particle index and n = 
1,2,- •■ ,27V is the Lyapunov index. Eqs. ©, © and 
(|13fl lead to the third property of the localization width 



yy(4JV-n+i) _ yy(")_ 



(14) 



This is the conjugate property for the localization width 
WW. 

Some of the advantages and disadvantages of the use 
of the localization width VV(*0 to discuss the Lyapunov 
localizations are as follows. An advantage of the localiza- 
tion width W( n ) is that its calculation does not require 
the thermodynamic limit, different from the localization 
length fi( n ) defined by Eq. ©. Beside, the localization 
width W( n ' is applicable to the Lyapunov localization 
seen in Fig. Hfb), whereas the localization length fiw 

requires that the tail of the quantity 7^ (t) decay expo- 
nentially in space. Moreover the localization width 
does not require an artificial parameter like the thresh- 
old 9 S (0, 1) in the quantity C^ n) defined by Eq. 0, 
so it can be determined from the dynamics only. On 
the other hand, a disadvantage of the localization width 
V\A n ) is that using this quantity, for example, we cannot 
distinguish between one peak of height a and width 2/3 
with constants a and /3 and two peaks of height a and 
width (3 which should be recognized as different localized 
behaviors. Therefore, in principle we have to look at the 
concrete shapes of localized peaks to check the relation 
between this localization width and the peak width in 
actual numerical calculations, even if we calculate the lo- 
calization width of the Lyapunov vector. The quantity 
C { q } defined by Eq. also has the same disadvantage. 



FIG. 3: Normalized localization widths W {n) /N in the case 
of d = 0.5. The circles, triangles and squares correspond to 
the cases of iV = 25, 50 and 100, respectively. The broken 
line interrupted by dots, the solid line and the broken lines 
(and the dotted line) correspond to W (n) = W wa v (~ 0.736iV), 
Wran(~ 0.6517V) and W m in(~ 2), respectively, (a) Full scale 
of the normalized localization widths as functions of the nor- 
malized Lyapunov index n/(2N). (b) An enlarged graph of 
the normalized localization widths corresponding to the small 
positive Lyapunov exponent region as functions of the Lya- 
punov index n. Symbols with a grey-filled background corre- 
spond to the two-point steps of the Lyapunov spectra which 
are indicated by similar circles in Fig. 2] Symbols surrounded 
by a rectangle of broken lines correspond to the four-point 
steps of the Lyapunov spectra which are indicated by similar 
rectangles surrounding symbols in Fig. [1] 



Fig. El are examples of graphs of localization widths 
as functions of Lyapunov index in quasi-one-dimcnsional 
systems. Here we have plotted the localization widths 
H , W /N normalized by the particle number N for quasi- 
one-dimensional systems of size TV = 25 (circles) , N — 50 
(triangles) and N = 100 (squares) where d = 0.5 (den- 
sity p = 0.523 • ■ • ). In order to take the time-average in 
Eq. (JSJ to calculate the localization width, we sample 

the quantities 7^ (t) In 7^ (t) just after collisions over 
1000A collisions and take arithmetic averages of them. 
(In this paper we always calculate localization widths 
in this way.) Fig. [Sfa) is the full scale graph of the 
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localization widths as functions of the normalized Lya- 
punov index n/(2N), and Fig. Etb) is a part of the 
localization widths corresponding to the small positive 
Lyapunov exponent region as functions of the Lyapunov 
index n. (Note that we use the normalized Lyapunov 
index n/(2N) for the x-axis in Fig. Efa), which is dif- 
ferent from the x-axis as the Lyapunov index n itself in 
Fig. [Jb).) To plot the values of the localization widths 
W (n \ j = 2N + 1,2N + 2, • • • , AN is omitted because of 
the conjugate relation 1)14(1. In the region of large pos- 
itive Lyapunov exponents the localization width W'™' 
is an monotonically increasing function of the Lyapunov 
index n. This implies that localized behavior of the Lya- 
punov vector is stronger in the large (in absolute value) 
Lyapunov exponent region. The shape of the localization 
width in this region is similar qualitatively to the shape 
given by the quantity defined by Eq. {7|) which 

was calculated in a square system [l^. Figure |3{ a) also 
shows that the value of the normalized Lyapunov local- 
ization widths y\A n J /N decrease as a function of particle 
number N in the region where Lyapunov spectra change 
smoothly. Noting that the localization width has a lower 
bound (W^ n '/N > 1/iV), this suggests the existence of 
the thermodynamic limit for the spectrum of the local- 
ization widths. 



FIG. 4: Lyapunov spectra normalized by the largest Lya- 
punov exponent A' 1 ^ as functions of the Lyapunov index n 
for the case of d = 0.5. The circles, triangles and squares 
correspond to the cases of JV = 25, 50 and 100, respectively. 
Symbols with a grey-filled background correspond to symbols 
with similar circles in the graph|3fb). Symbols surrounded by 
a rectangle of broken lines correspond to symbols surrounded 
by similar rectangles in the graph |3Jb). Inset: The full scale 
graph of the positive branch of the Lyapunov spectra as func- 
tions of the normalized Lyapunov index n/(2N). 



Another important point in Fig. |31 is the connection 
with the stepwise structure of Lyapunov spectra. Fig. 0] 
is the Lyapunov spectra normalized by the largest Lya- 
punov exponent A' 1 ' as a function of the Lyapunov index 
n in the case of N = 25 (circles) , N = 50 (triangles) and 
N = 100 (squares) with d = 0.5. These three cases corre- 
spond to the three cases of the localization widths in Fig. 
03 In this figure the stepwise structures of Lyapunov spec- 
tra appear in the small Lyapunov exponent region. The 
values of the largest Lyapunov exponents are w 1.28 
for N = 25, A^ « 1.31 for N = 50, A^ » 1.31 for 
N = 100. The inset of Fig. 0|is the full scale graphs of 
the positive branch of the Lyapunov spectra normalized 
by the largest Lyapunov exponents as functions of the 
normalized Lyapunov index n/(2N) (Note again that we 
use the different horizontal axes in the main graphs and 
the inset in Fig. 0J), and show that the global shapes of 
these graphs are similar and independent of the particle 
number N. Comparing Fig. 0] with Fig. |3{b) for the 



localization widths in the same region of the Lyapunov 
index n, we notice that the localization widths W'*'/JV 
normalized by the particle number N corresponding to 
the clear two-point steps of the Lyapunov spectra take 
almost same value which is almost independent of the 
sequence of the steps and the particle number N (See 
the symbols with a grey-filled background in Figs. 01b) 
and 0] about this point.). A similar thing can be seen 
for the four-point steps of the Lyapunov spectra. (See 
the symbols surrounded by a rectangle of broken lines in 
Figs. Elb) and^J) One may also notice that the region of 
the localization widths corresponding to the stepwise re- 
gion of the Lyapunov spectra, which is around the region 
n/(2N) > 0.8 of the Lyapunov index n in Figs. ap- 
proximately, can be distinguished from the other region. 
These points about the localization widths that connect 
with the stepwise structure of the Lyapunov spectra are 
other merits of the use of the localization width to discuss 
the behavior of the Lyapunov vectors. 



FIG. 5: Localized regions (grey-filled regions surrounded by 
solid lines) of Lyapunov vectors ^r' 1 ' and the numbers of col- 
liding particle pairs (as pairs of black-filled circles) as func- 
tions of the collision number n t and the particle number j in 
a quasi-one-dimensional system of N — 50 and d = 10° . The 
solid line is the contour plot 7^ = 0.3 of the amplitudes of 
the normalized Lyapunov vector components 7^ of the j'-th 
particle corresponding to the largest Lyapunov exponent A' 1 ', 
which corresponds to the three dimensional plot given in Fig. 

Eft). 



Before finishing this section it is valuable to summarize 
some specific values of the localization width W^™- 1 which 
have clear physical meanings. The first value of the lo- 
calization width is yy(") = Wmax = N which means 
that the amplitudes |<5r^™' ) | of the Lyapunov vector com- 
ponents for each particle take the same values as each 
other, as mentioned before. The Lyapunov localization 
close to this value actually occurs in one of the zero- 
Lyapunov exponent, as shown in W^^/N in Fig. [3fa). 
The second value is the lower bound of the Lyapunov lo- 
calization: VV ( ") > VVmin = 2. This is a little stronger 
condition than the inequality Ww > 1 which we have 
already shown in the inequality l|10|) . This value of the 
localization width is shown in Fig. [31 as the dotted line 
(the three lines of W {n) /N = 2/N = 0.08 for N = 25, 
2/N = 0.04 for N = 50, and 2/N = 0.02 for N = 100). 
This lower bound for the localization width comes from 
the fact that particle collisions occur between two parti- 
cles, and is partly supported by the fact that the width 

(n) 

of the the amplitudes 7^ of the normalized Lyapunov 
vector components of the j-th particle corresponding to 
the n-th Lyapunov exponent \^ is almost 2 in the low 
density limit and in large Lyapunov exponents, for exam- 
ple, as shown in Fig. Hfb). It should also be emphasized 
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that there is a relation between the localized region of a 
Lyapunov vector and positions of colliding particle pairs. 
To show this point, in Fig. [S] we plot colliding parti- 
cle pairs as well as a contour plot of the amplitude 
as functions of the collision number n t and the particle 
number j in a quasi-one-dimensional system of TV = 50 
and d — 10 5 , which correspond to the three dimensional 
plot of the Lyapunov localization given in Fig. |3b) . It is 
clear from this figure that particle collisions occur at the 
beginning of the sharp rectangular shapes of the localiza- 
(n) 

tion of the amplitudes 7^ . This suggests that the Lya- 
punov localization comes from a short interaction range 
of particle-particle interactions. The minimum value of 
the localization width will also be discussed in the sub- 
sections and [^21 in this paper. The third value of 
the localization width is the value for the case in which 
the amplitudes Jj (t), k — 1,2, ••■ ,N are distributed 
randomly with an equal probability (the solid line in Fig. 
EJ. It is given by M/W = W ran « 0.6517V approximately, 
and gives an upper bound on the localization widths of 
the Lyapunov index corresponding to the continuos part 
of the Lyapunov spectrum in many cases [2^ ■ The fourth 
value of the localization width is the value obtained for 
a wavelike structure in the Lyapunov vector. The bro- 
ken line in Fig. [3] is given as the localization width 
of 7j (n) (t) = ((|aj Tl) sin(27rjn/7V + (3)\ 2 )) where the con- 
stant is determined by the normalization condition 

of ^^{t) and ((■■•)) means the average over (3 randomly 
chosen by the flat distribution. Our numerical estimate 
says that the localization width W^ in such a case is n- 
independent and is given by W^ n ' = W wav ~ 0.736iV ap- 
proximately. Fig. Otb) show that some of the localization 
widths arc in this range (the broken line interrupted by 
dots) approximately, and they actually correspond to the 
transverse Lyapunov mode, namely a wavelike-structure 
of Lyapunov vectors corresponding to the two-point steps 
of the Lyapunov spectrum 0, 0] . 



III. DENSITY DEPENDENCE OF THE 
LYAPUNOV LOCALIZATION 



In this section we compare the Lyapunov localization 
widths Vy*-™-* as a function of the Lyapunov index at dif- 
ferent densities. We concentrate our consideration on the 
region where the Lyapunov spectra change smoothly as 
a function of the Lyapunov index, and show that there 
are two kinds of Lyapunov index dependence of the local- 
ization widths: exponential dependence [ST>] and linear 
dependence [CD] . In the remaining sections of this paper 
the particle number N — 50 is constant, and we change 
the particle density p by changing the parameter d con- 
nected to p by Eq. J2J) . 



A. Lyapunov Localization at High Density 



Figure is the Lyapunov localization widths W>' n ) /N 
normalized by the particle number N as functions of 
the normalized Lyapunov index n/(2N) in the cases of 
d = -0.05 (density p = 0.826 ••• , circles), d = 0.3 
(density p — 0.604- ••, triangles) and d = 1.5 (den- 
sity p = 0.314- ••, squares). In this figure the local- 
ization width W (n) = W ran ~ 0.657V in the case of ran- 
dom components and the minimum localization width 
yyX") = Wmin = 2 are indicated by the solid line and the 
broken line, respectively. 



FIG. 6: Normalized localization widths W' n ' /N as functions 
of the normalized Lyapunov index n/(2N) in high density 
cases of d — —0.05 (circles), d = 0.3 (triangles) and d = 1.5 
(squares). The data are fitted by exponential functions. The 
localization widths VV^ n ' = Wran and Wmin are indicated by 
the solid line and the broken line, respectively. The localiza- 
tion width indicated by the arrow corresponds to the localized 
behavior of the Lyapunov vector shown in Fig. ffla) . 



In Fig. [S] we fitted the graphs to exponential func- 
tions y = ad + fid exp^x) with fitting parameters 
otdi (3d and y^. Here we find the best values of the 
fitting parameters to be (a-o.05,/3-o.05, 7-0.05 ) = 
(0.627454,-0.366969,-13.1761), (a . 3 ,A).3,7o.3) = 
(0.590977,-0.469104,-9.33407) and (0:1.5,^1.5,71.5) = 
(0.572315,-0.476522,-5.70097). The localization 
widths yv^™) are nicely fitted by these exponential 
functions in the region n/(2N) < 0.6 (the exponential 
dependence [£T> ]), 

Figure also shows that the localization widths de- 
crease as the parameter d increases (therefore as the den- 
sity p decreases). The largest Lyapunov exponent is an 
increasing function of the particle density, as will be dis- 
cussed in Sect. |V] so this result means that the Lya- 
punov localization is stronger in a system with weaker 
chaos characterized by a smaller value of the largest- 
Lyapunov exponent. It should also be noted that in 
this region the localization widths are smaller than the 
localization width of random elements, which is about 
W r an ~ 0.65iV, and larger than W min = 0MN. Fig. 
Ufa) in Sect. UTI was the behavior of the Lyapunov local- 
ization corresponding to the localization width indicated 
by the arrow in Fig. |SJ It is important to note that 
the shape of the quantities Jj(t) as a function of the 
particle index j have a tail like that in Fig. Ufa) in the 
region showing an exponential dependence [£ T>\ of the 
localization widths. 
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B. Lyapunov Localization at Low Density 

Now we consider the Lyapunov localization in low den- 
sity cases, given by large values of the parameter d. 



FIG. 7: Normalized localization widths W' n ' /N as functions 
of the normalized Lyapunov index n/(2N) in low density cases 
of d = 10 2 (squares), d = 5 x 10 2 (triangles) and d = 10 5 
(circles). The localization widths W' n ' = W ra n and Wmin 
are indicated by the solid line and the broken line, respec- 
tively, (a) Localization widths in the Lyapunov index region 
n/ (2N) < 0.7. The data are fitted by an exponential function. 
The grey region is the region in which the localization widths 
show a linear dependence [CD] as a function of Lyapunov in- 
dex in the low density limit, (b) Enlarged graphs including 
a linear dependence region of the localization widths. The 
data are fitted by linear functions. This figure also includes a 
part of the graph of the normalized localization widths in the 
case of d = 1.5 (the crosses), which has already been shown in 
Fig. HJ The localization widths indicated by the black arrow, 
the grey arrow and the outline arrow correspond to localized 
behaviors of the Lyapunov vectors shown in Figs. Pfa). Pfb) 
and|H| respectively. The thick grey horizontal line to connect 
the localization widths indicated by the black arrow and the 
outline arrow in Fig. (b) is given to show that these values 
take almost the same value. 



this paper. (In Fig. |7f a) the region of the linear de- 
pendence [CD] in the low density limit is shaded grey.) 
Secondly, in Fig. Ufa), we can still recognize a region of 
localization widths in which the Lyapunov index depen- 
dence of the Lyapunov localization widths is exponential 
(the dependence [££>])• To show it clearly we give a fit 
of an exponential function y — a' + (3' exp(7'x) with val- 
ues (a', j3',i) = (0.560814,-2.73348,-7.50543) of the 
fitting parameters in Fig. |7f a). It should be emphasized 
that the shapes of localization widths are almost inde- 
pendent of the particle density p, at least in the three 
density cases shown in Fig. Ufa). As the third point 
one may notice that from above characteristics of the lo- 
calization widths in the low density cases the region of 
the linear dependence [CD] and the region of the expo- 
nential dependence [£D] are more distinguishable in a 
lower density case, with an accompanying sharp bending 
of the localization width profile. Figure [7\ suggests that 
the spectrum of the localization widths as a function of 
Lyapunov index has a distinct shape in the low density 
limit p — > with a critical value (el of the Lyapunov 
index where the localization width shows a linear 

dependence [CD] in the Lyapunov index n < (el and 
shows an exponential dependence [£D] in the Lyapunov 
index n > (el- 



Figure |7{a) is the Lyapunov localization widths 
yy(™) /TV normalized by the particle number N as func- 
tions of normalized Lyapunov index nj (2N) in the case of 
d = 10 2 (density p = 0.00777 • • • , squares), d = 5 x 10 2 
(density p — 0.00156 • • ■ , triangles) and d = 10 5 (den- 
sity p = 0.00000785 ••• , circles), which correspond to 
the region where the Lyapunov spectra change smoothly. 
It is clear that the Lyapunov index dependences of the 
localization widths in this figure are different from the 
high density cases shown in Fig. E|in a couple of senses. 
First, linear dependences of the localization widths with 
respect to Lyapunov index appear in the region of small 
localization widths, which correspond to large Lyapunov 
exponents. In Fig. EJt>), fits of Lyapunov index de- 
pendences of the localization widths by linear functions 
y = &dX+/3d are given with fitting parameters ay and (3d'- 
(<5 10 2,/3 102 ) = (0.503874, 0.046405) in the case of d= 10 2 , 
(5 5x i 2,/3 5x i 2) = (0.347437,0.044468) in the case of 
d = 5 X 10 2 , and (a loB ,/3 10 5) = (0.224129,0.0427506) 
in the case of d = 10 5 . It is important to note that 
the localization widths are always larger than 2, namely 
W {n) /N > 2/N = 0.04 in this figure, and the smallest lo- 
calization widths corresponding to the largest Lyapunov 
exponents are close to this minimum value in these low 
density cases. Besides, graphs of the linear dependence 
of the localization widths is flatter in the lower density 
case, as shown by the fact that the value of the fitting 
parameter ad decreases as the value of the parameter d 
increases. The existence of the linear dependence [CD] 
of the localization widths is one of the main results of 



(15) 

FIG. 8: Localized behaviors of the amplitudes 7] of the 
normalized Lyapunov vector components of each particle as 
a function of the collision number nt and the particle index 
j in the case of d = 5 x 10 2 . The corresponding localization 
width is indicated by the outline arrow in Fig. |7|Lh 



It is important to note that the exponential depen- 
dence [ED] and the linear dependence [CD] of the local- 
ization widths can also be distinguished by the shapes of 
the amplitudes 7^ (t) of the normalized Lyapunov vector 
components of each particle. To discuss this point it is 
important to note that Fig. 12 a) is the localized behavior 
of the Lyapunov vector corresponding to the Lyapunov 
width indicated by the black arrow in Figs. Eland[7Jb) 
and is in the region of exponential dependence [£ D] of 
localization widths. This localization shows a long tail 
behavior. On the other hand Fig. [2Jb) is a Lyapunov 
localization corresponding to the localization width indi- 
cated by the grey arrow in Fig. H^b) and is in a region of 
linear dependence [CD] of localization widths. This lo- 
calization shows a sharp rectangular shape (with width 
2). These observations lead to the conjecture that in the 
exponential dependence [£D\ the amplitudes jj (t) have 
a long tail, and in the linear dependence [CD] the am- 
plitudes Jj^fy) have a rather sharp rectangular shape. 
To make this point convincing, we compare the behav- 
iors of the amplitudes 7^ (t) corresponding to two local- 
ization widths which take almost the same value of the 
localization widths but show different dependences [CD] 
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and [£T>] of the localization widths. For this purpose we 
choose the 15-th localization width W^ 15 ) for d = 5 x 10 2 . 
This localization width, indicated by the outlined arrow 
in Fig. Efb), is in the region of the linear dependence 
[CD] and takes almost the same value as the localization 
width W^ 1 ) indicated by the black arrow in Figs. EJand 
Efb) in the case d — 1.5, which is in the region of expo- 
nential dependence [£T>] of the localization widths, and 
whose localized behavior is shown in Fig. Ufa). Figure |H1 

is the localized behavior of the amplitudes Jj (t) as a 
function of the collision number n t and the particle index 
j in the case d — 5 x 10 2 . Here, the data in Fig. [5] is 
taken after every second collision. This figure shows that 
some peaks of the amplitudes jj (t) have fiat tops with 
a width 2 and distinct rectangular shapes, although its 
corresponding Lyapunov exponent, the 15-th Lyapunov 
exponent is far from the largest one and it is less localized 
than in Fig. 

It should be noted that the linear dependence [CD] of 
the localization widths as a function of Lyapunov index is 
not a property of the quasi-one-dimensional systems only. 
In Appendix IB 21 we consider the localization widths in 
a square system, and show that the linear linear depen- 
dence [CD] of the Lyapunov localization also appears in 
the square system. 



IV. LYAPUNOV SPECTRA, AND THE ANGLES 
AND AMPLITUDES OF LYAPUNOV VECTOR 
COMPONENTS 

In this section, we consider how differences between 
the linear [CD] and exponential dependences [£T>] of the 
localization widths affect other quantities, such as the 
Lyapunov spectrum, the angle between the spatial and 
momentum parts of the Lyapunov vectors, and the am- 
plitudes of the x and y-components of the spatial part 
and the momentum part of the Lyapunov vectors. We 
also investigate the effects of the stepwise structure of 
the Lyapunov spectra on these quantities. 



exponents. The Lyapunov index dependence of the local- 
ization widths in these cases were shown in Figs. E3 and 
[7| It should be emphasized that in the low density cases, 
in which the linear dependence [CD] appears, we can rec- 
ognize a sharp bend in the Lyapunov spectra around the 
normalized Lyapunov index n — Cel- Such a bending 
point corresponds to the point that distinguishes the de- 
pendences [CD] and [£ V] of the localization widths in 
Fig- U\ (In this figure the region of the linear depen- 
dence [CD] of localization widths at low densities has a 
grey background.) It may be noted that such a bend- 
ing of the Lyapunov spectrum has also been reported in 
a fully two-dimensional system |2!| (Also see the Lya- 
punov spectra for the Fermi- Pasta- Uhmi models in Refs. 



FIG. 9: Lyapunov spectra normalized by the largest Lya- 
punov exponents A' 1 ' as functions of normalized Lyapunov 
index n/(2JV) in the cases of d = —5 x 10~ 2 (white-filled cir- 
cles), d — 0.3 (white-filled triangles), d = 1.5 (white-filled 
squares), d = 10 2 (black-filled squares), d — 5 x 10 2 (black- 
filled triangles) and d = 10 5 (black-filled circles). The grey 
region is the region in which the localization widths show a 
linear dependence [CD] as a function of Lyapunov index at 
low densities. The small figure in the upper right side is the 
enlarged graphs of the small positive region of the Lyapunov 
spectra for the low density cases of d= 10 2 , 5 x 10 2 and 10 5 . 



In Fig. El we can see some stepwise structures of the 
Lyapunov spectra, not only in the high density cases d = 
—5 x 10 -2 , d = 0.3 and 1.5, but also in the low density 
cases d = 10 2 , d = 5 x 10 2 and d = 10 5 . (See the enlarged 
graphs in the upper right side of Fig. |5]for the stepwise 
structure of the Lyapunov spectra in these low density 
cases.) The steps of the Lyapunov spectra consist of two- 
point steps and four point steps, which were considered in 
Ref. [2l| in detail. Different from the high density cases, 
at low density the separations of the two-point steps and 
the four-point steps in the Lyapunov spectra are not so 
clear. 



A. Lyapunov Spectra 

Figure is the Lyapunov spectra normalized by their 
largest Lyapunov exponents AW for the systems as func- 
tions of the normalized Lyapunov index n/(2N) for d — 
-5 x 1CT 2 (white-filled circles), d = 0.3 (white-filled tri- 
angles), d — 1.5 (white-filled squares), d — 10 2 (black- 
filled squares), d = 5 x 10 2 (black-filled triangles) and 
d = 10 5 (black-filled circles). Here the values of the 
largest Lyapunov exponents are given by A' 1 -* ~ 4.79 for 
d = -5 x 10~ 2 , AW « 1.64 for d = 0.3, A« w 0.769 
for d = 1.5, AW « 0.0579 for d = 10 2 , AW « 0.0161 for 
d = 5 x 10 2 , A (1) « 0.000157 for d = 10 5 . The negative 
branch of the Lyapunov spectra are omitted in this figure, 
because of the conjugate property (|12|) of the Lyapunov 



B. Angle Between the Spatial and the Momentum 
Parts of the Lyapunov Vectors 

As a second example to illustrate the effects of the 
two dependences [CD] and [£T>] of localization widths, 
we consider the angle 8^ between the spatial part 
5q(") = (£qj , ? ' ' ' i^Ijv) 7 an d th e momentum 
part Sp( n ^ — (Sp^ , 6p% , ■ ■ ■ , 5pff) T of the Lyapunov 
vector 6T^ n ' corresponding to the n-th Lyapunov expo- 
nent \( n \ which is defined by 



l ) = cos" 1 



SqM . g p (r. 
|£q(")| |<5p( n )| 



(15) 
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Noting the conjugate relations 

j q (4JV-n+l) , ^ p (4jV-n+l) ^ j q (n) , fi p (n) 
|j q (4JV-n+l)| |,5p(4Ar-ri+l)| ~ _ |£ q (n)| |<5p(™)| ^ ' 

following from the conjugate property (| 1 1 ft of the Lya- 
punov vectors, the angle flW satisfies the condition 

e (4iV-n+l) = 7r ( 17 ) 

This angle has already been investigated in many-particle 
systems [3, |2l| . It should be noted that some analytical 
approaches give formula to calculate Lyapunov exponents 
from the spatial part only (or the momentum part only) 
of the tangent vector [jj, 132, H3 , so it is important to 
investigate the relation between the spatial part and the 
momentum part of the Lyapunov vector. 



FIG. 10: Time-average (6^ n ')/Tv of the angle divided by 7r 
for the spatial part <5q(™' and the momentum part 5p' n ' of 
the Lyapunov vector corresponding to the n-th Lyapunov ex- 
ponent A*™- 1 as functions of the normalized Lyapunov index 
n/{2N) for d = -5 x 1CT 2 (white-filled circles), d = 0.3 
(white-filled triangles), d = 1.5 (white-filled squares), d = 10 2 
(black-filled squares) , d = 5 x 10 2 (black-filled triangles) and 
d = 10 (black-filled circles). The grey region is the region in 
which the localization widths show a linear dependence [CD] 
as a function of Lyapunov index in low density cases. The 
solid line corresponds to the value 8^ — ty/2 of the angle. 
The small figure in the upper right side is the enlarged graphs 
for the low density cases d — 10 2 , 5 x 10 2 and 10 5 in the region 
of large Lyapunov indices. 



In Fig. we show the graphs of the time-average 
(6>( n ))/7r of the angle 0W divided by 7T as functions 
of the normalized Lyapunov index n/(2N) in the cases 
d = — 5 x 10~ 2 (white-filled circles), d = 0.3 (white- 
filled triangles), d = 1.5 (white-filled squares), d = 10 2 
(black-filled squares), d = 5 x 10 2 (black-filled trian- 
gles) and d = 10 5 (black-filled circles). Here the time- 
average of the angle 9^ is the arithmetic average of 
their values immediately after collision, for lOOOiV col- 
lisions. The plots of the time-averaged angles (9^), 
n = 2N + 1, 27V + 2, • • • , 47V, corresponding to negative 
branch of the Lyapunov exponents are omitted because 
of the conjugate property (|T7Jl of the angle 6^ n K It is 
noted that in our numerical calculations the amplitudes 
|(Sp (n) |, n = 27V - 2, 27V - 1 and 2N corresponding to 
zero-Lyapunov exponents are zero (or extremely small), 
so the angles 9 (n \ n = 2N - 2, 2N - 1 and 27V cannot 
be defined. This is the reason why we do not plot these 
angles in Fig. (See the subsection IIV CI about this 
point.). This figure shows that in the low density cases, 
in which the linear dependence [CD] of the localization 
widths appears, the spectra of the averaged angle (9^) 



bends around the value n/(2N) = Cbl/(27V) of the nor- 
malized Lyapunov index. (In Fig. HOI we showed a region 
of the linear dependence [CD] of the localization widths 
as a grey region.) In the region of the linear dependence 
[CD] of the localization widths the time-averaged angles 
increase monotonically as functions of Lyapunov 
indices n, while the time-averaged angles (9^) corre- 
sponding to the exponential dependence [£D] of the lo- 
calization widths is close to tt/2, meaning that in this 
region the vector <5q( n ) is almost orthogonal to the vec- 
tor tfpOO (See the case of d= 10 5 in Fig. HOI ). 

We can also see the effect of the stepwise structure of 
the Lyapunov spectra in the angle 9^ n \ As shown in Fig. 
^| the values of the time-averaged angles (9^) corre- 
sponding to the two-point steps of the Lyapunov spectra 
in Fig. El are rather smaller than the values correspond- 
ing to their four-point steps, in all densities shown in Fig. 
El (See the upper right side of Fig. ^Jfor the low density 
cases). Therefore the Lyapunov index dependence of the 
time-averaged angles (6^ n >) can be used to distinguish the 
two-point step of the Lyapunov spectra from the other 
parts. 



C. Amplitudes of the x and y-Components of the 
Lyapunov Vectors 

As shown in the subsection IIV Al in low density cases 
the positive branch of Lyapunov spectra is separated into 
three parts: (i) A region in which the Lyapunov spectrum 
is a rapidly decreasing function of Lyapunov index and 
corresponds to the linear dependence [CD] of the local- 
ization widths, (ii) A region in which values of the Lya- 
punov exponents are very small compared to the largest 
Lyapunov exponent and corresponds to the exponential 
dependence [£D] of the localization widths, and (hi) A 
region showing the stepwise structure of the Lyapunov 
spectra. The boundary of the regions (i) and (ii) of the 
Lyapunov spectra appear as a sharp bending of the Lya- 
punov spectra in the low density limit. To understand 
this characteristic of the Lyapunov spectra we note that a 
positive Lyapunov exponent is connected to a decorrela- 
tion rate in chaotic dynamics, and a larger positive Lya- 
punov exponent means that a trajectory diverges more 
rapidly on the energy surface in ergodic systems and loses 
correlation with its initial value more quickly. This is 
partly supported by the fact that the Kolmogorov-Sinai 
entropy, which is equal to the sum of all the positive Lya- 
punov exponents (Pesin's identity) in closed systems, is 
connected to the inverse relaxation time in ergodic sys- 
tems 0. A relation between a Lyapunov exponent and 
a decay rate of time-correlation is also discussed in Refs. 
|34l l35| . Another example of a relation between a time- 
scale of the dynamics and a positive Lyapunov exponent 
is the tracer particle effect in Lyapunov spectra [36| . On 
this subject a many-particle system consisting of a small 
light particle called the tracer particle moving at high 
speed (the short time scale) and many big heavy particles 
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moving at low speeds (the long time scale) is considered, 
and it is shown that the existence of the tracer particle 
leads to relatively large Lyapunov exponents compared 
with the Lyapunov exponents corresponding to the other 
heavy particles. The above considerations suggest that 
at low density the quasi-one-dimensional systems exhibit 
a bending of the Lyapunov spectra, which separates the 
set of positive Lyapunov exponents into two groups (i) 
and (ii) of relatively different values, and this suggests 
the existence of a time-scale separation in the dynamics 
of these systems. 

One might think that in the quasi-one-dimensional sys- 
tems such a time-scale separation may come from the 
very narrow rectangular shape of the system. In the 
quasi-one-dimensional system, especially in the large d 
case with a low density, the time-scale of oscillation of 
each particle in the x-direction is much larger than the 
time-scale of that in the y-direction, and it may be one 
possibility leading to the time-scale separation mentioned 
above. (However this possibility may be rather unlikely, 
because as shown in Appendix IB II the bending point of 
the Lyapunov spectra does not change significantly by 
changing the system from the quasi-one-dimcnsional sys- 
tem to the square system with the same area L x L y and 
the same number of particles N.) Another possibility 
to explain this bending of the Lyapunov spectrum may 
be from the different roles of the spatial and momentum 
parts of the Lyapunov vectors at low density. The result 
in the previous subsection supports this point, showing 
that the spatial and momentum parts of the Lyapunov 
vectors are not in the same direction at low density. 

Motivated by the above considerations, in this subsec- 
tion we consider the amplitudes of the four parts of the 
Lyapunov vectors: the x and y-components in the spatial 
part and the momentum part of the Lyapunov vectors, 
namely the four kinds of quantities defined by 



k=l 



«5T< n) (*) 



(18) 



(19) 



j — x and y, in each Lyapunov index n. Here Sq^} (t) 
and 5p$(t) are the j-component of the spatial coordi- 



nate part Sq^ 



= {SQxfc {t), SQyj. 1 (t)) and the momentum 



part 5p£' = (5p^fe (t),5p^ ; (*)) of the Lyapunov vector 
_ (£q(™) ( Sp^ ) corresponding to the fc-th particle 
and the n-th Lyapunov exponent at time t, respectively. 
These four quantities are normalized as 



by their definitions. Using the conjugate property i|ll|) 
of the Lyapunov vector we obtain 



(AN-n+l) 



nfr n+l) (t)=r> { :ht). 



(21) 
(22) 



Therefore we can omit the quantities T]^ N n+1 ^ (t) and 



(20) 



^(4JV n+i)^ n _ 2, ■ • • , 2N corresponding to the 
negative branch of the Lyapunov spectra. 



FIG. 11: Time- averages of the normalized amplitudes of the 
re-component of the spatial part ((r/qx (t)) , circles), the y- 
component of the spatial coordinate part ( (rj^y 1 (t)), triangles) , 
the ^-component of the momentum part ({rjpx(t)}, crosses), 
and its y-component of the momentum part ({i]py\t)}, di- 
agonal crosses) of the Lyapunov vectors as functions of the 
normalized Lyapunov index n/(2N) for d = —5 x 10 -2 (the 
graph (a)), d — 0.1 (the graph (b)), d = 1 (the graph (c)) 
and d = 10 5 (the graph (d)). The small panels in the top 
right of figures (a), (b), (c) and (d) are enlarged graphs of 
(?7pS (t)} and {f]py {t)} in small positive Lyapunov exponent 
regions. The grey region in the graph (d) is the region in 
which the localization widths show a linear dependence [CD] 
as a function of Lyapunov index. 



Figure^Jis the graphs of the time-averaged quantities 
(Vqx\t)) (circles), (rj$ '(*)) (triangles), (?7^(f)) (crosses) 
and (r)py\t)) (diagonal crosses) of the quantities defined 
by Eqs. l|18|l and 119|l as functions of the normalized 
Lyapunov index n/(2N) for d = —5 x 1CT 2 (Fig. ITTT a). 
density p = 0.826 •••), d = 0.1 (Fig. HiTbl. density 
p = 0.713 ■■■), d=l (Fig. Etc), density p = 0.392 • • • ) 
and d = 10 5 (Fig. UTTdL density p = 0.00000785 ••• )■ 
Here the time-average {X) of the quantities X = r)qx (t) , 
r]qy\t) r]p$ (t) and rjpy(t) is taken as the arithmetic av- 
erage of the quantity X immediately after collisions, for 
1000N collisions. The cases of d = -5 x 10~ 2 and d = 10 5 
shown in Figs. lUf a! and II If d) are the cases of the high- 
est density and the lowest density considered in Sect. IIIII 
and the previous two subsections of this section, respec- 
tively, and the other two cases are given to show interme- 
diate situations between the two cases of Figs. Illf a) and 
lllf dh It is noted that Figs. Illf ah (b) and (c) are for the 
case of a particle density in which the linear dependences 
[CD] of localization widths do not appear yet, and only in 
the grey region of Figs. Illf d) does the linear dependence 
[CD] appear. From this figure it is clear that there is 
a strong asymmetry between the spatial and momentum 
parts of the Lyapunov vectors in any density region, as 
well as an asymmetry between the x and y-components 
of the Lyapunov vectors especially at low density. 

The behavior of the graphs in Fig. ^2 are different in 
the region where the Lyapunov spectra change smoothly 
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and in the stepwise region of the Lyapunov spectra. 
First we discuss the region where the Lyapunov spectra 
change smoothly, which is the region of the Lyapunov 
indexes n/(2N) < 0.8 approximately in Fig. ^] In the 
high density case shown in Fig. Illf a) the momentum 

parts ({rjpx (t)) and (%}(£))) of the Lyapunov vectors 
are larger than the corresponding spatial parts ({rjqx (i)) 

and (r/qy^ (t))). As the density decreases, the spatial parts 
of the Lyapunov vectors increase (and momentum parts 
decrease), first in the large positive Lyapunov exponent 
region (See Fig. Illf al — > (b).), and then in the small Lya- 
punov exponent region (See Fig. Illf bl — ■> (c) — » (d).). 
At very low density, as shown in Fig. Illf dL the spa- 
tial parts of the Lyapunov vectors are much larger than 
the corresponding momentum parts for Lyapunov vec- 
tors of any index, and a linear dependence [CD] of the 
localization widths appears. In this density region a sig- 
nificant gap in the x and y-components of the spatial part 
of the Lyapunov vectors appears, and the transverse part 
({Vqy (t)) ) is larger than the longitudinal part (rjq 7 ^ (t)} in 
the large positive Lyapunov exponent region while while 
the reverse is true in the small positive Lyapunov expo- 
nent region. In the low density limit the momentum parts 
(Vpx\t)) and (r]py\t)} are extremely small so the spatial 
parts (r]qx (t)) and (rj q y (t)) are almost symmetric in the 
line y = 0.5 as required by the the normalization condi- 
tion (|2U|I . It should be emphasized that in the Lyapunov 
index of the linear dependence [CD] of the localization 

widths the quantities (%/ (*)) and (r) q ™ (t)) are rather 
flat, shown as the grey region in Fig. Illf d). We can show 
that these asymmetries between the x and y-components 
of the Lyapunov vectors actually come from its narrow 
rectangular shape, because as shown in Appendix IB 41 
such asymmetries do not appear in the square system, 
although the asymmetries between the spatial and mo- 
mentum parts of the Lyapunov vectors still exist there. 

As a second point, Fig. ^2 also shows some char- 
acteristics of the time-averaged amplitudes of the Lya- 
punov vector components corresponding to the stepwise 
region of the Lyapunov spectra. In Fig. ^2 we can see 
flat parts consisting of two-points (four-points) for the 
time-averaged normalized amplitudes (T]qx{t)}, (Vqy\t)}i 
(r]px CO) and (rjpy\t)) of the Lyapunov vectors, which 
are corresponding to the two-point steps (the four-point 
steps) of the Lyapunov spectra (See Figs. ITTT al. (b), 
(c) and (d) in a region of the Lyapunov index around 
n/(2N) > 0.9, as well as the small panels in the top right 
of these figures for enlarged momentum parts (rjpx (t)) 
and (ripy (t)) in the small positive Lyapunov exponent 
regions.). Values of the two-point steps of the transverse 
part (rj q y (t)) and (rjpy (i)) of the Lyapunov vectors cor- 
responding to the two-point steps of the Lyapunov spec- 
tra are larger than their longitudinal parts (r]qx(t)) and 
(VpxHt)) in Figs. ITTTaL (b), (c) and (d). On the other 
hand, values of the four-point steps of the longitudinal 



part (i]qx (t)) and {ry P x (t)) of the Lyapunov vectors cor- 
responding to the four-point steps of the Lyapunov spec- 
tra are larger than their transverse parts (rjqy\t)) and 

(Vpy(i)) m Figs. HTT a - ). (b) and (c), and the spatial parts 
of the Lyapunov vectors in Figs. Illf dL At any density 

the spatial parts (%" (£)) and (%,(*)) are rather larger 
than the momentum parts (f]px(t)) and (r]py\t))i respec- 
tively, in the region of Lyapunov indices indicating the 
stepwise structure of the Lyapunov spectra. It should be 
emphasized that wavelike structure (the transverse Lya- 
punov modes) in the y-components of the spatial part 
and the momentum part of the Lyapunov vectors ap- 
pears in the two-point steps of the Lyapunov spectra re- 
ported by and Fig. ITT1 suggests that we may get 
a rather clear wave-like structure in the spatial part of 
the Lyapunov vectors than in its momentum part. We 
can conclude a similar result for the four-point steps of 
the Lyapunov spectra in which time-depending wavelike- 
structure in the y-components of the spatial part of the 
Lyapunov vectors and in cc-c omp onents of the spatial 
part of the Lyapunov vectors |15l l2lj . It may be noted 

that the time-averaged amplitudes (r/p 7 ^ (t)} and (rjpy (i)), 
n = 2N — 2, 2N — 1 and 2N corresponding to zero- 
Lyapunov exponents are almost zero at any density (See 
the small panels in the top right of Figs. ITTT aL (b), (c) 
and (d).). 



V. DENSITY DEPENDENCE OF THE 
LARGEST LYAPUNOV EXPONENT AND 
RELATED QUANTITIES 

As we have already shown in Sect. IIIII the linear de- 
pendence [CD] of localization widths as a function of Lya- 
punov index appears in cases of low density, but it was 
not so clear at what density the linear dependence [CD] 
of the Lyapunov localization starts to appear as the par- 
ticle density p deceases from 1 . In this section we discuss 
this problem by calculating density dependences of the 
largest Lyapunov exponent, the angle and amplitudes of 
Lyapunov vector components, and the localization width, 
which correspond to the largest Lyapunov exponent. Es- 
pecially, we show that the particle density region in which 
the linear dependence [CD] appears is almost the same 
as the density region in which the largest Lyapunov ex- 
ponent begins to satisfies the Krylov relation. 



A. Largest Lyapunov Exponent 

Fig. El is the largest Lyapunov exponent as a 
function of the density p in the quasi-one-dimensional 
system. In this figure we showed the numerical error 
e = A (1) + A< 4Ar ) (which must be zero in the Hamiltonian 
system) at each data point as an error bar of the length 
2 1 e | , showing that the values of the largest Lyapunov ex- 
ponents at very low densities are less accurate than those 
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at higher densities. The grey region in this figure is the 
density region in which the linear dependence [CD] of the 
localization widths appears. 



ntlf, and assume that the mean free path If is inversely 
proportional to the density p, so we obtain 



FIG. 12: Largest Lyapunov exponent A' 1 ' as a function of 
particle density p. The grey region is the density region in 
which the linear dependence [CD] of the localization widths 
as a function of Lyapunov index appears. Numerical data 
is fitted by the Krylov relation A' 1 ' ~ ap In (f3p) (solid line) 
and the power function y = a' p 13 (broken line) with fitting 
parameters a, j3, a 1 and 0' . The error bar in each data point 
is given by the absolute value 2|A' X ' + A' 4iY '| of twice of the 
sum of the largest and smallest Lyapunov exponent. 



A' 1 



-u th jp ]n(Ryp). 



(25) 



using a constant 7 of If ~ (7/5) 1 . Comparing Eq. H25|) 
with Eq. J23J) we obtain the relation 



a 



Uth 

R ' 



(26) 



It is known that in the low density limit the largest 
Lyapunov exponent A*- 1 -* should have the form 



A (1) ~ ap ln(/3p) 



(23) 



with parameters a and (3. Eq. i|23|) is called the Krylov 
relation j3?J and has already been demonstrated numer- 
ically in a fully two-dimensional system consisting of 
many-hard-disks H^H^, apart from a factor. In Fig. ^| 
we fitted the numerical data to the Krylov relation i|23|) 
with parameter values a = —1.66875 and (3 — 1.28252 
(the solid line), which gives a good fit for the density 
dependence of the largest Lyapunov exponent. It is im- 
portant to note that the density region that satisfies the 
Krylov relation almost coincides with the density region 
in which the linear dependence [CD] of the localization 
widths appears. To test the Krylov relation l|23|l we 
also tried to fit the numerical data by a power function 
y = a' (broken line) with fitting parameters a' and (3' 
in Fig. We used the parameter values a 1 = 2.55126 
and f3' = 0.808158, and it is shown in Fig. U^that the 
Krylov relation (|23|l gives a better fit than this power 
function in the grey region. 

Now we check the values of the fitting parameters a 
and (3 used to fit the graph of Fig. ^| by following the 
rough derivation of the Krylov relation in Ref . 1] . (Note 
that more exact derivations are known for some specific 
systems, such as the Lorentz gas, etc. See, for example, 
Refs. |39|,|40|-) First we note that after nt particle col- 
lisions the amplitude of a Lyapunov vector is stretched 
by a factor (lf/R) n * approximately with If the mean free 
path. Introducing the collision rate v = n t /t to connect 
the mean free path If with the time t we estimate the 
largest Lyapunov exponent as 



lim - In ( — 

t \ R 



= vm—. 
R 



(24) 



We approximate the collision rate v by a thermal velocity 
u t h = \J (2/M)(K/N) as v ~ Uth/lf because of u t ht ~ 



which is independent of the value of the constant 7. 
Using the parameter values u t h = V% and R = 1 in 
our numerical calculations, we obtain a/f3 ~ — v2 = 
— 1.4142 • • • , which is consistent of the values of the fit- 
ting parameters a and (3 in Fig. a/0 = —1.3011 



FIG. 13: Mean free time 77 as a function of density p. The 
dotted line is a plot of the fitting function y = 1/ (u t h'yx) with 
the thermal velocity u t h and the fitting parameter 7. The grey 
region is the density region in which the linear dependence 
[CD] of the localization widths as a function of Lyapunov 
index appears. 



Next we proceed to estimate the factor 7 = (Ifp) _1 . 
Fig. 1131 is the graph of the mean free time Tf as a func- 
tion of the density p. Noting the relation If ~ u t hJth, in 
this figure we fitted the numerical data by the function 
y = 1/ (uth'yx) with the value 7 = 13.341 for the fitting 
parameter 7. (It may be emphasized that this function 
gives a nice fit in the density region in which the linear 
dependence [CD] of the Lyapunov localization appears.) 
However this value 7 = 13.341, which is (3 in Eq. Ij23(l in 
the case of R = 1, is about 10 times as large as the fit- 
ting value of /3 used to fit the graph of Fig. It should 
be noted that because of the quasi-one-dimensional prop- 
erty of the system, a well known relation If ~ l/(\/2p) 
[41| for fully two-dimensional hard-disk systems with a 
Maxellian velocity distribution cannot give a correct re- 
lation in the quasi-one-dimensional system, although this 
gives a parameter value (3 which is rather close to the 
value used in Fig. It may be noted that in order 

to derive the expression l|25() for the largest Lyapunov 
exponent we neglected some characteristics of the quasi- 
one-dimensional systems, such as the fact that in the 
quasi-one-dimensional system with a large value of the 
parameter d particles collide mainly head on and it is 
rather rare for particles to have grazing collisions with 
other particles. We should take these points into account 
to get more precise expressions for the parameters a and 
13 in Eq. (J23J than in Eq. (H>J). 
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B. Angle and Amplitudes of Lyapunov Vector 
Components 

As the next example we consider the angle 9^ defined 
by Eq. I|15|l between the spatial and momentum parts of 
the Lyapunov vector corresponding to the largest Lya- 
punov exponent. Figure 1141 is the graph of the time- 
average of the angle 9^ /n divided by 7r as a function 
of particle density p. The grey region in this figure is 
the density region in which the linear dependence [CD] 
of the localization widths appears. The graph has a lo- 
cal maximum point about p ~ 0.2, which is close to the 
value of particle density in which the linear dependence 
[CD] of the localization widths starts to appear as the 
density decreases from 1. As we have already discussed 
in Sects. [H] and II I II the existence of the linear depen- 
dence [CD] of the localization widths can be checked by 
a linear dependence of the localization widths as a func- 
tion of Lyapunov index and a rectangular shape of the 
amplitude of the Lyapunov vector, but it is rather hard 
to use this criterion to distinguish the density region of 
the linear dependence [CD], because it is not easy to 
recognize these behaviors in an intermediate region be- 
tween the density region of the exponential dependence 
[£D] only and the density region of both the dependences 
[ED] and [CD]. In this sense the graph of the angle 
as a function of particle density p may give a more dis- 
tinct criterion to distinguish the density region in which 
the linear dependence [CD] appears. It is noted that the 
angle 6^' does not seems to go to zero in the low density 
limit as shown in the inset of Fig. ITTI 



is dominant and a gap appears between its x-component 
{viqx{t)) and the y-component (r) q ])(t)). On the other 
hand, in a high density region in which the linear de- 
pendence [CD] does not appear, the momentum parts 
(f]px (t)) and {ripy '(£)) of the Lyapunov vector are domi- 
nant. It should be noted that in the low density in which 
the linear dependence [CD] appears, the transverse com- 
ponents (riqy(t)) and (r)py (t)) are always larger than the 

longitudinal components {rfqx (t)) and (?jpi?(i)), respec- 
tively. 



FIG. 15: Time-averages of the normalized amplitudes of 
the i-component of the spatial part {(rjqx(t)}, circles), the 
y-component of the spatial part ((r]qy (t)}, triangles), the x- 
component of the momentum part ((r)px (i)), crosses), and 
the y-component of the momentum part ((rj^y (t)), diagonal 
crosses) of the Lyapunov vector corresponding to the largest 
Lyapunov exponent as functions of the particle density p. The 
grey region is the density region in which the linear depen- 
dence [CD] of the localization widths as a function of Lya- 
punov index appears. 



C. Localization Width 



FIG. 14: Density dependence of the time-average (9^}/tt 
of the angle divided by n between the spatial part and the 
momentum part of the Lyapunov vector corresponding to the 
largest Lyapunov exponent as a function of the particle den- 
sity p. The grey region is the density region in which the 
linear dependence [CD] of the localization widths as a func- 
tion of Lyapunov index appears. The inset: The same graph 
but including the angle 6^ at lower density. 



Next we consider the density dependence of the time- 
average of normalized amplitude of the x-component and 
the y-component in the spatial part and the momentum 
part of the Lyapunov vector corresponding to the largest 
Lyapunov exponent, which are defined by Eqs. (|18J) and 

(fH^I in n = 1. Fig. El is the graphs of (rjqx \t)} (cir- 
cles), {rjqy(t)) (triangles), (rjpj(t)) (crosses) and (rj^(t)) 
(diagonal crosses) as functions of the particle density p. 
Again the density region in which the linear dependence 
[CD] of the localization widths appears is indicated as a 
grey region. This may not be a good criterion to distin- 
guish the density region of the linear dependence [CD] , 
but it is clear that in the density region of the linear de- 
pendence [CD] the spatial parts of the Lyapunov vector 



As the last quantity in this section, we consider the 
density dependence of the Lyapunov localization width 
corresponding to the largest Lyapunov exponent. Fig- 
ure 1161 is the graph of the Lyapunov localization width 
W" /N normalized by the particle number TV as a func- 
tion of density p. This figure shows that in the low den- 
sity limit the normalized localization width W^ 1 - 1 /N goes 
to the value 2/N(= 0.04) (solid line in Fig. |TBJ), which 
is discussed as the minimum localization width W m in/N 
in Sect. [HI In Fig. ^|we indicate the density region of 
the linear dependence [CD] as the grey region, but it is 
not so clear from this figure how to distinguish the den- 
sity region in which the linear dependence [CD] of the 
localization widths appears. 



FIG. 16: Normalized localization width W' 1 ' /N correspond- 
ing to the largest Lyapunov exponent as a function of density 
p. The solid line is the minimum value of the normalized lo- 
calization width: Wmin/N = 0.04. The grey region is the 
density region in which the linear dependence [CD] of the 
localization widths as a function of Lyapunov index appears. 
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VI. CONCLUSION AND REMARKS 

In this paper we have discussed localized behaviors of 
Lyapunov vectors (the Lyapunov localization) for quasi- 
one-dimensional systems consisting of many hard disks 
(with periodic boundary conditions except in Appendix 
IB 5|) . The quasi-one-dimensional system was introduced 
as a particle system whose shape is a very narrow rectan- 
gle that does not allow the exchange of particle positions. 
We compared some methods to characterize the localized 
behavior of the Lyapunov vectors, and as one of such 
methods in this paper we used a quantity called the lo- 
calization width, whose logarithm is given by an entropy 
for the amplitude distribution of the Lyapunov vector 
components of each particle. The localization width indi- 
cates the number of particles contributing to the localized 
part of Lyapunov vector components. It could not only 
be used as an indicator to measure the magnitude of the 
localized behavior of the Lyapunov vectors, but also be 
used to distinguish different delocalized properties of the 
Lyapunov vectors such as the derealization associated 
with a random distribution of particle component ampli- 
tudes, a derealization associated with a uniform distri- 
bution and a derealization associated with a wave-like 
structure (corresponding to stepwise structure of the Lya- 
punov spectra ) in the Lyapunov vector (Figs. l^and^J. 
The localized region of the Lyapunov vectors is related to 
the points of colliding particles (Fig. |SJ), and this leads to 
the lower bound 2 of the localization width (Figs. [7|and 
I16p . Using the localization width we showed that there 
are two kind of the Lyapunov localizations in many-hard- 
disk systems. The first type of the Lyapunov localization 
is characterized by an exponential dependence [£T>] of 
the localization widths as a function of Lyapunov index 
(Fig. ISJ , and by its long tale of localized Lyapunov vec- 
tors (Fig. Efa))- This type of the Lyapunov localization 
is observed at any particle density. The second type of 
the Lyapunov localization is characterized by the linear 
dependence [CD] of the localization widths as a function 
of Lyapunov index (Fig. UJ, and by the sharp rectangu- 
lar shape of the localized Lyapunov vectors (Figs. Hfb) 
and|SJ). This type of Lyapunov localization appears only 
in low density cases and in Lyapunov indices correspond- 
ing to the large Lyapunov exponents in absolute values. 
We showed that in the density region of the linear de- 
pendence [CD] of the localization widths the Lyapunov 
spectra is bent and separated into two parts (except for 
the stepwise region of the Lyapunov spectra): one cor- 
responding to the exponential dependence [£D] and take 
very small values corresponding to the largest Lyapunov 
exponent, and the other corresponding to the linear de- 
pendence [CD] and shows a rapid decreasing dependence 
of Lyapunov index (Fig. [§J|. It was also shown that differ- 
ences between the exponential dependence [ED] and the 
linear dependence [CD] appear in the angle 9^ between 
the spatial and momentum parts of the Lyapunov vec- 
tors (Fig. I1U[I and in the amplitudes of the ^-component 
and y-component of the spatial part (Fig. IHf d)). (Here 



we took the y-direction as the narrow direction of the 
rectangle and the x-direction as the longer orthogonal 
direction.) The density region, in which the linear depen- 
dence [CD] of localization widths appears, almost exactly 
coincides with the density region in which the density de- 
pendence of the largest Lyapunov exponent A^ 1 ^ satisfies 
the Krylov relation (Fig. H2)l . We also indicated that 
at the boundary of the density region of the exponential 
dependence [£T>] only and the density region of both the 
linear dependence [CD] and the exponential dependence 
[£T>], the angle 9^ corresponding to the largest Lya- 
punov exponent shows a local maximum as a function of 
particle density (Fig. 114(1 . 

In this paper we observed differences in the amplitudes 
of the x and y-components of the Lyapunov vectors (Figs. 
ITT1 and 115(1 . These differences come from the difference 
in the roles of the directions in the quasi-one-dimensional 
systems. We also observed differences in the amplitudes 
of the spatial and momentum parts of the Lyapunov 
vectors. In the region where the Lyapunov spectra are 
changing smoothly the amplitude of the spatial part of 
the Lyapunov vectors is larger (smaller) than of the mo- 
mentum part in low (high) density cases (Figs. 1111 1151 
and I2U(I . The spatial and momentum parts of the Lya- 
punov vectors are in almost the same direction at high 
density, whereas they are rather close to orthogonal in 
the low density case, especially in the region of the expo- 
nential dependence [£T>] of the localization widths as a 
function of Lyapunov index (Figs. H0land lT5|) . These be- 
haviors are found not only in the quasi-one-dimensional 
systems but also in the square system (Fig. |3UJ|. Con- 
cerning the stepwise region of the Lyapunov spectra, as 
shown in Fig. 1111 the amplitudes of the y-component, 
(transverse component), of the spatial and momentum 
parts of the Lyapunov vectors are larger than the corre- 
sponding ^-components, namely (their longitudinal com- 
ponents) in the two-point steps of the Lyapunov spectra, 
whereas they are opposite in four-point steps of the Lya- 
punov spectra (except in very low density cases like in 

pig. una)). 

The particle number dependence of the Lyapunov lo- 
calization width normalized by the particle number N 
was investigated a little in this paper. Figure [3] suggests 
that the Lyapunov localization width normalized by the 
particle number decreases as a function of the particle 
number itself, but we need a further calculation to show 
the existence of its thermodynamical limit [N — » +oo), 
as well as its low density limit (p — ► 0). 

We suggested that the bending of the Lyapunov spec- 
tra in the low density cases accompanying the linear de- 
pendence [CD] of the localization widths may come from 
a time-scale separation of the dynamics. As one of the 
possible to explanations of this point we demonstrate a 
strong asymmetry of the amplitudes of the spatial and 
momentum parts of the Lyapunov vectors at low den- 
sity. However this may not give an enough explanation 
as to where the bending point of the Lyapunov spectra 
is, and one may also indicate that the Lyapunov spec- 
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trum is related to growing (or reducing) speeds of the 
amplitudes of the Lyapunov vectors, not the amplitudes 
of the Lyapunov vectors themselves. This remains an 
open question. 

Results in this paper suggest that the localized behav- 
ior of the Lyapunov vectors comes from the short range 
property of particle interaction (Fig. ©, noting that the 
hard-core interaction of the systems used in this paper 
is the shortest range interactions possible. Therefore it 
should be interesting to compare our results with the 
Lyapunov localization of systems with long range parti- 
cle interactions. Moreover, if the Lyapunov localization is 
connected to the existence of the thermodynamical limit 
of the largest Lyapunov exponent [l4|, then it should 
be interesting to investigate the Lyapunov localization in 
the Fermi-Pasta-Ulam model in which a numerical cal- 
culation suggests a logarithmic divergence of the largest 
Lyapunov exponent as a function of particle number N 
|42| . These remain as future problems on the subject of 
Lyapunov localizations. 
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T(x) = xhxx — x + 1. 



(A2) 



It is easy to show that the function T(x) satisfies the 
inequality 



T{x) > 0, in x > 0. (A3) 

Eqs. © and (|A1|) . and the inequality (|A3|I lead to 



In TV - S {7 



N' 



N 



7 <">(t)2V 



> 0, (A4) 



3=1 



noting the inequality 7 j"^(t)iV > 0. Therefore we ob- 
tain the inequality W^™-* < N using Eq. ©. From the 
two inequalities 1 < W (n) and W {n) < N we derive the 
inequality l|10l) for the localization width. 

We can also show that the function T(x) defined by Eq. 
(IA2I) satisfies T(x) — only when x = 1. Using this point 
and Eqs. © and (|A4|) we get the fact that the equality 
yy(«) — ]\[ for the localization width is satisfied only 
when all the quantities 7 j (£), j = 1, 2, • • • N are equal, 



namely when 7^ (t) — j [ - z> (t) 



„( 2 ) / 



1/N. 



APPENDIX B: COMPARISON WITH SQUARE 
CASES 



APPENDIX A: INEQUALITY FOR THE 
LOCALIZATION WIDTH 



In this appendix we derive the inequality (|10(l for the 
localization width of the Lyapunov vector. 

First we note the inequality < which comes from 
Eqs. © and ©. This leads to the inequality 1 < W (n) . 

Second, we note 



N 



3=1 



= N~ 



AT 



N 



In 



3=1 



In 



3=1 



(Al) 



where we used Eq. @, and the function T{x) of x is 
defined by 



In this appendix we discuss differences between quasi- 
one-dimensional systems and square systems for the Lya- 
punov localization and its related phenomena. For the 
quasi-one-dimcnsional system in this paper we chose the 
system lengths as 



L y = L' y = 2R(l + 10~ 6 ), (Bl) 
L x = L' x = NL y (l + d) (B2) 
[in the quasi-one-dimensional cases] 

with a parameter d to change the particle density. For 
meaningful comparisons between the square and quasi- 
one-dimensional cases we use system lengths for the 
square systems so that both cases give the same area, 
for the same value of the parameter d, namely 



L, 



L'U 



= v/JV(l + d)[2fl(l + 10- 6 )] 2 
[in the square cases]. 



(B3) 



Except for these lengths L x and L y we use the same 
parameter values as given in the text, such as R = 1, 
M = 1 and E = N . In this appendix we consider systems 
of 50 particles (N = 50). In Appendixes lBl1lB2llFT3l and 
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IB 41 we consider the case of d = 10 5 in which the linear 
dependence [CD] of localization widths with respect to 
Lyapunov index appears. In Appendix IB 51 we discuss 
boundary effects in the localization width in the case of 
d = 0.5. 



1. Lyapunov Spectra 

The first example is the Lyapunov spectra of the quasi- 
one-dimensional case and the square case. Fig. El is 
the positive branch of the Lyapunov spectra normalized 
by the largest Lyapunov exponent A^ 1 ' as functions of 
the normalized Lyapunov index n/(2N) in the quasi-one- 
dimensional system (circles) and the square system (tri- 
angles) in the case of d — 10 5 . Here the values of the 
largest Lyapunov exponents are given by AW w 0.000157 
in the quasi-one-dimensional system, and A^ « 0.000552 
in the square system. As shown in Fig. 1171 a sharp bend- 
ing of the Lyapunov spectra occurs in both cases, but it 
occurs in a slightly smaller Lyapunov index in the square 
case than in the quasi-one-dimensional case. The Lya- 
punov spectrum in the quasi-one-dimensional case is also 
shown in Fig. El 



FIG. 17: Lyapunov spectra normalized by the largest Lya- 
punov exponent A^ 1 ' as functions of the normalized Lya- 
punov index n/(2N) in the case of d — 10 in the quasi- 
one-dimensional system (circles) and the square system (tri- 
angles). Inset: Enlarged graphs of the normalized Lyapunov 
spectra in a region of small positive Lyapunov exponents. 



There is a difference between the quasi-one- 
dimensional case and the square cases in the region 
of small positive Lyapunov exponents. In the inset 
of Fig. El we can recognize two stepwise structures 
of the Lyapunov spectrum consisting of six Lyapunov 
exponents in the quasi-one-dimensional system. It 
should be noted that each of these steplike structures 
consists of one four-point step and one two-point step 
although these two steps are too close to be distinguished 
in Fig. El (This is also evidenced by investigating 
wavelike structures of the Lyapunov vectors, the so 
called Lyapunov modes |2l|.) On the other hand we 
cannot recognize such a stepwise structure in the square 
case, because 50 particles in the square system with 
the density p = 0.00000785 • • • are too small to show 
stepwise structure of the Lyapunov spectrum. (It is 
known that rectangular systems including the quasi-one- 
dimensional systems show a longer stepwise region of 
Lyapunov spectra than in the square system with the 
same area pj.) In the square system the gap between the 
smallest non-zero positive Lyapunov exponent A2jv-4 
and the zero Lyapunov exponents is larger than in the 
quasi-one-dimensional system. 



2. Localization Widths 

Next we consider the localization width in the quasi- 
one-dimensional case and the square case. Fig. 1181 is 
the localization widths W^/N, n = 1, 2, • • • , 2N nor- 
malized by the particle number N as functions of the 
normalized Lyapunov index n/(2N) in the quasi-one- 
dimensional case (circles) and the square case (triangles) 
in the case of d = 10 5 . In the region of Lyapunov 
index corresponding to the Lyapunov spectra chang- 
ing smoothly, qualitative behavior of the localization 
width are rather similar, and consist of the linear de- 
pendence [CD] in small Lyapunov indices and the ex- 
ponential dependence [£D] in other indices, although 
the region of the linear dependence [CD] in the square 
case is slightly smaller than in the quasi-one-dimensional 
case. In Fig. I18j _we give only fits of exponential func- 
tions y — a + /3exp{7x} for the exponential depen- 
dence [£D] by dotted lines. (The parameter values to 
fit data are (a, /3, 7) = (0.560814,-2.73348,-7.50544) 
for the quasi-one-dimensional case, and (a, $, 7) = 
(0.597656, -3.35386,-10.9674) for the square case.) The 
localization widths are larger than their minimum value 
yyw = Wrnin = 2 and are smaller than its random com- 
ponent case W (n) = W ran ~ 0.651N except in the step- 
wise region of the Lyapunov spectra, and the smallest 
localization width are close to W m in in both the 

cases. 



FIG. 18: Normalized localization width W <n) /N as func- 
tions of the normalized Lyapunov index n/(2N) in the case 
of d = 10 5 in the quasi-one-dimensional case (circles) and 
the square case (triangles). The dotted lines are fits by ex- 
ponential functions (the dependence [£2?]). The broken line 
interrupted by dots, the solid line and the broken line cor- 
respond to W (n) = W wav (fst 0.736AO, Wran(~ 0.651.2V) and 
W m m(= 2), respectively. 



As in the Lyapunov spectra, the localization width 
shows a different behavior with the Lyapunov index cor- 
responding to the stepwise structure of the Lyapunov 
spectra. As we have already mentioned, in the square 
case of N — 50 there is no stepwise region of the Lya- 
punov spectrum, so there is not a corresponding structure 
in the localization width. We can see that in the square 
case the localization widths W^ n ' decrease as a function 
of the Lyapunov index n in the region of large Lyapunov 
indices corresponding to small positive Lyapunov expo- 
nents. On the other hand we can see the effect of the step- 
wise structure of the Lyapunov spectrum on the localiza- 
tion width of the quasi-one-dimensional case. Especially 
a pair of two dots on the line W'"' = W wav correspond 
to the two-point step of the Lyapunov spectrum shown 
in Fig. 1171 supporting that in the Lyapunov spectrum 
of Fig. El f° r the quasi-one-dimensional case the first 
(second) step of the Lyapunov spectrum is a four-point 
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step (a two-point step) looking from the zero-Lyapunov 
exponents. 

3. Angle Between the Spatial and Momentum 
Parts of the Lyapunov Vector 

Fig. O is the time-averaged angles {9^)/w, n = 
1, 2, • • • , 2N divided by tt between the spatial and mo- 
mentum parts of the Lyapunov vectors in the cases of 
the quasi-one-dimensional system (circles) and the square 
system (triangles) as functions of the normalized Lya- 
punov index n/(2N) for d = 10 5 . In both the cases, the 
angle 0^ n ' is a rapidly increasing function of Lyapunov 
index n in the linear dependence [CD] region of the lo- 
calization width, and almost tt/2 in the exponential de- 
pendence [£D] region. 



FIG. 19: Time-average (9^)/ty of the angles divided by 7r 
between the spatial and momentum parts of the Lyapunov 
vectors in the quasi-one-dimensional case (circles) and the 
square case (triangles) as functions of the normalized Lya- 
punov index n/(2N) for d — 10 5 . The solid line corresponds 
to the value 6^ n > — tt/2 of the angle. Inset: Enlarged graphs of 
normalized angles (9^)/% in the region of Lyapunov index 
in which the the exponential dependence [£T>] of the local- 
ization widths (and the stepwise structure of the Lyapunov 
spectrum) appears. 



We can see the stepwise structure in the angle 6^ n ' 
as a function of Lyapunov index n in the quasi-one- 
dimensional case (See the inset of Fig. E3)- This struc- 
ture corresponds to the stepwise structure of the Lya- 
punov spectrum shown in Fig. 1171 



4. Amplitudes of the x- and {/-Components of the 
Spatial and Momentum Parts of the Lyapunov 
Vectors 

As the next quantities to investigate the difference 
between the quasi-one-dimensional case and the square 
case, we consider the amplitudes of the x- and y- 
components of the spatial and momentum parts of the 
Lyapunov vectors. Fig. |20 is the graphs of the time- 
averages of normalized amplitudes of the a;-component of 
the spatial part ((rjq X (t)), circles), the y-component of 
the spatial part {(r)qy\t)), triangles), the x-component 
of the momentum part ((f]^; (t)}, crosses), and its 
y-component ({f]^y(t)}, diagonal crosses) of the Lya- 
punov vectors as functions of normalized Lyapunov in- 
dex n/(2N) in the square case for d — 10 5 . The corre- 
sponding graphs for the quasi-one-dimcnsional case are 
shown in Fig. Illf dh In the square system there is no 
difference between the x- and y-directions, so there is no 



gap between the x- and y-components in the spatial part 
and the momentum part of the Lyapunov vectors (except 
for those where the spatial part corresponds to the zero- 
Lyapunov exponents). Although the momentum parts 

(Vpx(t)) an d (Vpv (t)) (f ne spatial parts and 
(riqy (£))) are rapidly decreasing (increasing) functions of 
Lyapunov index in the region of the linear dependence 
[CD] of the localization widths (See the small figure in 
the top right of Fig. [2U for the momentum parts of the 
Lyapunov vectors. Also note the normalization condi- 
tion (|2(Jfl to know the corresponding spatial parts.), the 
boundary of the two dependences [CD] and [£T>] is not 
so clear in this figure. 



FIG. 20: Time- average of the normalized amplitudes of the 
a;-component of the spatial part ((7?g2 (£)), circles), the y- 
component of the spatial part ((rjqy\t)}, triangles), the x- 
component of the momentum part ((ffc™ (t)), crosses), and 
the y-component of the momentum part ((?7py (i)) , diagonal 
crosses) of the Lyapunov vectors as functions of normalized 
Lyapunov index n/(2N) in the square case for d — 10 5 . Cor- 
responding quasi-one-dimensional case is given in Fig. Illf dh 
The small panels in the top right are the graphs of (r/p"\t)) 
and (^(t)) only. 



In Fig. 1201 we cannot recognize a stepwise structure of 

the quantities (r) q x (t)), (%} (<)}, fajS (*)) and (Vpy{t)}, 
because there is no stepwise structure of the Lyapunov 
spectrum in this system. It may be noted that the 

normalized x-component (r]p$ (t)} and the normalized y- 

component (rjpy (t)} of the momentum part of the Lya- 
punov vectors corresponding to the zero-Lyapunov expo- 
nents \( 2N ~ 2 \ xi 21 ^- 1 ) and \( 2N *> are almost zero, like in 
the quasi-one-dimensional case. 

5. Boundary Effects in the Localization Widths 

As the last example in this appendix we check bound- 
ary effects in Lyapunov localizations. In this subsection 
we calculate the localization widths for square systems 
and quasi-one-dimensional systems with purely hard-wall 
boundary conditions or purely periodic boundary con- 
ditions. To give reasonable comparisons between the 
the hard-wall boundary cases and the periodic boundary 
cases we choose the system lengths L x and L y so that 
the effective region for a particle to move is the same in 
both cases. This means that if we choose the system size 
as L x = L x and L y = L y for the periodic boundary case 
then they must be L x = L x + 2R and L y = L y + 2R for 
the corresponding hard-wall boundary case. Except for 
this point, we chose the same value of the other parame- 
ters for these two boundary cases. In this subsection we 
use the value d = 0.5 for the parameter d. Note that in 
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this case the linear dependence [CD] of the localization 
widths does not appear. 



FIG. 21: Normalized localization widths /N as fac- 
tions of the normalized Lyapunov index n/(2N) for d — 0.5 
in a square system with periodic boundary conditions (cir- 
cles), in a square system with hard- wall boundary conditions 
(triangles), in a quasi-one-dimensional system with periodic 
boundary conditions (crosses), in a quasi-one-dimensional sys- 
tem with hard- wall boundary conditions (diagonal crosses). 
The dotted lines are fits of the numerical data by exponential 
functions. The broken line interrupted by dots, the solid line 
and the broken line correspond to W {n) = Wwav(& 0.736iV), 
Wran{~ 0.65LZV) and W mi „(= 2), respectively. 



Figure PHI is the localization widths W^/N, n = 
1,2, ••■ , 2N normalized by the particle number N in 
the case of the square system with the periodic bound- 
ary conditions (circles), the square system with the 
hard- wall boundary conditions (triangles), the quasi-one- 
dimensional system with the periodic boundary condi- 
tions (crosses), the quasi-one-dimensional system with 
the hard- wall boundary conditions (diagonal crosses). 
Fig. 1211 shows that the localization widths in the hard- 
wall boundary conditions are smaller than in the periodic 
boundary conditions in both of the square cases and the 
quasi-one-dimensional cases. However it should be noted 
that differences in the values of the localization widths 
are rather small in the region of small Lyapunov indices 
n in which the Lyapunov vectors are strongly localized. 
This is a natural result because the Lyapunov localization 
is a local behavior so that it should not depend strongly 
on global conditions such as the boundary conditions, ft 
may be noted that in very narrow systems like the quasi- 



one-dimensional system a strong boundary effect appears 
because particles can collide with walls much more often 
than in square cases, and this could be the reason why 
the differences in the values of the localization widths 
from the different boundary conditions are larger in the 
quasi-one-dimensional system than in the square system. 
A part of the localization widths as a function of Lya- 
punov index are nicely fitted by the exponential function 
y = a + (3 cxp{7x} for the exponential dependence [ED] 
with fitting parameters a, /3 and 7, shown as the dotted 
lines in Fig. |21 Here, the values of fitting parameters 
are chosen as (a, £,7) = (0.653629,-0.377495,-10.9742 
for the square case with the periodic boundary con- 
ditions, (a,/3, 7 ) = (0.638696,-0.422529,-10.3554) for 
the square case with the hard-wall boundary condi- 
tions, (a,/3, 7 ) = (0.578304,-0.464099,-8.41439) for 
the quasi-one-dimensional periodic boundary case, and 
(a, j9, 7) = (0.534472,-0.449403,-6.76248) for the 
quasi-one-dimensional hard-wall boundary case. It may 
also be noted that the part of the localization widths cor- 
responding to the Lyapunov spectrum changing smoothly 
seems to be slightly larger than the value W ra „ for 
the localization width for random components in the 
square system with the periodic boundary conditions and 
d = 0.5 in Fig. EU 

In the quasi-one-dimensional periodic boundary case in 
Fig. |^we can recognize two pairs of localization widths 
on the line VV^ = W wav corresponding to the two-point 
steps of the Lyapunov spectra, but there is not a localiza- 
tion width of such a value in the hard- wall boundary case. 
This comes from the fact that in the periodic boundary 
case the system satisfies total momentum conservation 
leading to the two-point steps of the Lyapunov spectrum 
whereas the system with hard-wall boundary conditions 
does not satisfy such a conservation in any direction (De- 
tail of this point was discussed in Ref. |2ll|.). 
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